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Mammalian WTAP is a regulatory subunit of the RNA 
N6-methyladenosine methyltransferase 
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The methyltransferase like 3 (METTL3)-containing methyltransferase complex catalyzes the N6-methyladenosine 
(m6A) formation, a novel epitranscriptomic marker; however, the nature of this complex remains largely unknown. 
Here we report two new components of the human m6A methyltransferase complex, Wilms' tumor 1 -associating 
protein (WTAP) and methyltransferase like 14 (METTL14). WTAP interacts with METTL3 and METTL14, and is 
required for their localization into nuclear speckles enriched with pre-mRNA processing factors and for catalytic ac- 
tivity of the m6A methyltransferase in vivo. The majority of RNAs bound by WTAP and METTL3 in vivo represent 
mRNAs containing the consensus m6A motif. In the absence of WTAP, the RNA-binding capability of METTL3 is 
strongly reduced, suggesting that WTAP may function to regulate recruitment of the m6A methyltransferase complex 
to mRNA targets. Furthermore, transcriptomic analyses in combination with photoactivatable-ribonucleoside-en- 
hanced crosslinking and immunoprecipitation (PAR-CLIP) illustrate that WTAP and METTL3 regulate expression 
and alternative splicing of genes involved in transcription and RNA processing. Morpholino-mediated knockdown 
targeting WTAP and/or METTL3 in zebraflsh embryos caused tissue differentiation defects and increased apoptosis. 
These findings provide strong evidence that WTAP may function as a regulatory subunit in the m6A methyltransfer- 
ase complex and play a critical role in epitranscriptomic regulation of RNA metabolism. 
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Introduction 

Nitrogen or oxygen atoms in RNA are ubiqui- 
tously methylated with S-adenosylmethionine (SAM, 
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Adomet) serving as the methyl donor [1-3]. The N6- 
methyladenosine (m6A) modification is one of the most 
frequent methylations in eukaryotic mRNA, accounting 
for over 80% of all RNA base methylations and, it has 
been observed in various species [4-8]. Recently, two 
independent studies combining m6A immunoprecipita- 
tion with high-throughput analysis revealed that the m6A 
modification tends to mainly occur in intragenic regions. 
The frequency is particularly high in the 3 '-end of cod- 
ing sequence (CDS) and the first quarter of the 3'-UTR, 
with an average of 1 m6A per 2 000 ribonucleotides. The 
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m6A modification occurs in highly conserved regions 
with the consensus sequence, RRACH (R = G or A; H = 
A, C or U) [9-13]. During brain development a dynamic 
change in m6A levels on RNA has been observed [10], 
suggesting the existence of RNA m6A demethylase(s). 
Indeed, two AlkB dioxygenase family proteins, fat mass 
and obesity associated gene (FTO) and ALKBH5, were 
recently demonstrated to catalyze the removal of the 
m6A mark both in vitro and in vivo [14, 15]. These two 
enzymes function in various biological processes such 
as development, metabolism, and fertility. They regulate 
the expression levels of thousands of genes, suggesting a 
pivotal epitranscriptomic function of m6A modification 
in regulating RNA metabolism [8, 16-19]. 

In contrast to the discovery of RNA m6A demethyl- 
ases, the nature of the methyltransferase complex respon- 
sible for catalyzing m6A formation in RNA still remains 
largely unknown. The mammalian m6A methyltransfer- 
ase complex is generally perceived to contain at least two 
key subcomplexes designated MT-A (200 kD) and MT-B 
(800 kD) [20, 21]. So far only a 70-kD protein, methyl- 
transferase like 3 (METTL3) has been identified as one 
component of the 200-kD MT-A complex. METTL3 ho- 
mologs have been identified in S. cerevisiae, Drosophila 
and Arabidopsis [22-26]. Although mRNA methylation 
occurs primarily within the RRACH consensus sequence 
in mammals, only a portion of the RRACH sites contain 
the m6A modification [9, 10, 27]. Mechanistic under- 
standing of the selectivity, regulation and function of the 
m6A modification depends on a complete characteriza- 
tion of the methyltransferase complex. 

In this study, we identify for the first time, Wilms' 
tumor 1 (WTl)-associating protein (WTAP) and meth- 
yltransferase like 14 (METTL14) as components of the 
mammalian m6A methyltransferase complex. In particu- 
lar, WTAP appears to serve as a regulatory subunit essen- 
tial for m6A methyltransferase activity. 

Results 

Identification of WTAP and METTL14 as METTLS-asso- 
ciating proteins 

To identify mammalian METTL3 -associating pro- 
teins, we utilized tandem affinity purification to pull 
down stably expressed SFB-tagged METTL3 from 293T 
cells. Subsequent mass spectrometry analysis of affinity- 
purified SFB-METTL3 protein complexes identified 
two potential METTL3 -associating proteins, WTAP and 
METTL14 (Supplementary information, Figure S1A). 
Consistently, AtFIP37 and Mum2, the WTAP homologs 
in plant and yeast, respectively, have previously been 
shown to associate with METTL3 in these species [26, 



28]. METTL14 protein shares 43% identity (Supple- 
mentary information, Figure SIB) with METTL3, and 
a previous phylogenetic analysis has suggested that the 
two proteins are homologs [23]. Furthermore, METTL14 
contains similar motifs essential for catalytic activity as 
METTL3 (Supplementary information, Figure SIB). The 
interactions among METTL3, WTAP and METT14 were 
confirmed by co-immunoprecipitation (Figure 1A, IE 
and Supplementary information, Figure S2A and S2B). 
The specificity of the interaction between METTL3 and 
WTAP was verified by examination of the endogenous 
interaction between these two proteins after depletion 
of WTAP by small interfering RNAs (siRNA) (Figure 
IB). GST pull-down using purified recombinant GST- 
METTL3 and 6><His-WTAP proteins demonstrated a di- 
rect interaction between METTL3 and WTAP (Figure 1 C 
and Supplementary information, Figure SIC). By using 
WTAP constructs with either the C- (WTAP-N, 1-200 
aa) or N-terminus (WTAP-C, 201-396 aa) truncated, we 
found that the interaction between WTAP and METTL3 
or METTL14 depends on the N-terminus of WTAP (Fig- 
ure ID and IF). It has been reported that RNA is not 
an essential component of the m6A methyltransferase 
complex [20, 29]. Consistently, we observed that nei- 
ther RNase A treatment nor inhibition of methylation by 
cycleoleucine affected the interaction between WTAP 
and METTL3 (Supplementary information, Figure S2C- 
S2E). Based on these findings, we conclude that WTAP 
and METTL14 are bona fide components of the m6A 
methyltransferase complex, and that both RNA and the 
m6A modification are dispensable for the interaction 
between WTAP and METTL3. In the remaining of this 
paper, we will refer to this complex as the WMM (WTAP, 
METTL3 and METTL14) complex. 

WTAP is required for m6A methyltransferase activity in 
vivo 

WTAP does not harbor any obvious catalytic domains, 
and in contrast to METTL3 and METTL14 that were 
recently demonstrated to form an active complex capable 
of catalyzing m6A formation, WTAP showed no catalytic 
activity itself or effect on the activity of METTL3-MET- 
TL14 complex in vitro [30]. Given the direct interaction 
between WTAP and METTL3, we speculated that WTAP 
may serve to regulate m6A methyltransferase activity 
in vivo. To test this hypothesis, we examined the m6A 
levels in mRNA extracted from WTAP- or METTL3- 
knockdown cells. Dot-blot assays showed that depletion 
of either of these proteins resulted in significantly lower 
m6A levels in mRNA (Figure 2A and Supplementary in- 
formation, Figure S3A). RT-PCR and western blot analy- 
ses verified the knockdown efficiency and specificity, 
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Figure 1 WTAP interacts with METTL3 and METTL14. (A) 293T cells were transfected with Flag-WTAP and Myc-METTL3 
constructs as indicated. Forty-eight hours later, cells were lysed and the lysates were subjected to immunoprecipitation using 
anti-Myc (Myc-IP) followed by immunoblotting with the anti-Flag antibodies. (B) 293T cells were treated with control siRNA 
(siCTRL) or siRNA targeting WTAP (siWTAP) for 48 h. Then cells were lysed and the lysates were subjected to IP using anti- 
WTAP. The immunoprecipitated samples were analyzed by immunoblotting with the anti-METTL3 antibodies. (C) Purified 
recombinant His-WTAP proteins were mixed with either GST or GST-METTL3 proteins as indicated, pulled down with GST 
beads, and subjected to immunoblotting with the indicated antibodies. (D) 293T cells were co-transfected with Myc-METTL3 
and Flag-WTAP full-length (-FL), N-terminal (-N) or C-terminal (-C) constructs as indicated. Forty-eight hours later, cells were 
lysed and the lysates were subjected to Myc-IP followed by immunoblotting with the anti-Flag antibodies. (E) 293T cells were 
transfected with Flag-WTAP and HA-METTL14 constructs as indicated. Forty-eight hours later cells were lysed and the ly- 
sates were subjected to HA-IP followed by immunoblotting with the anti-Flag antibodies. (F) 293T cells were co-transfected 
with HA-METTL14 and Flag-WTAP full-length (-FL), N-terminal (-N) or C-terminal (-C) constructs as indicated. Forty-eight 
hours later, cells were lysed and the lysates were subjected to HA-IP followed by immunoblotting with the anti-Flag antibod- 
ies. Supportive data were included in Supplementary information, Figures S1 and S2. 
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respectively (Supplementary information, Figure S3B 
and S3C). Moreover, depletion of WTAP or METTL3 
did not affect the expression of the other complex com- 
ponents (Supplementary information, Figure S3B and 
S3C). WTAP was originally identified as a protein asso- 
ciated with WT1 [31]. We examined whether WT1 might 
also be required for m6A methyltransferase activity in 
vivo. Interestingly, we found that WT1 depletion had no 
influence on m6A abundance, suggesting that WTAP's 



role in regulating m6A modification is independent of 
its previously described binding partner WT1 (Figure 
2A and Supplementary information, Figure S3A). Taken 
together, these findings suggest that WTAP affects m6A 
methyltransferase activity in vivo. 

WTAP regulates accumulation of METTL3 and 
METTL14 in nuclear speckles 

Previous studies have indicated that m6A modification 
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Figure 2 WTAP regulates the nuclear speckle localization of METTL3 and METTL14. (A) The graph represents the quantifi- 
cation of three independent dot-blot experiments (raw data were included in Supplementary information, Figure S3A). The y- 
axis represents the relative intensity of dots relative to that of the control group. P values were calculated using a two-tailed t- 
test. Error bars represent SD. (B-D) After transfection (48 h) with the indicated fluorescence FAM-labeled siRNA, HeLa cells 
were fixed and immunostained with the indicated antibodies. DNA was visualized by DAPI. Scale bar, 10 ^m. Supportive data 
were included in Supplementary information, Figure S3. 
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occurs at the pre-mRNA stage within the nucleus and it 
thereby may affect splicing and mRNA export [2, 32- 
35]. Nuclear speckles are enriched in proteins involved 
in RNA processing and alternative splicing. Interest- 
ingly, METTL3 and WTAP as well as the m6A demeth- 
ylases FTO and ALKBH5 have previously been shown 
to colocalize with nuclear speckle markers [14, 15, 31]. 
The physical interaction between WTAP and METTL3 
and METTL14 prompted us to examine the cellular 
distribution of METLL14. Using immunofluorescence 
microscopy, we found that endogenous METTL14, like 
METTL3 and WTAP, exhibits a diffused nucleoplasmic 
staining pattern with intense granule-like foci, and that 
it colocalizes well with the pre-mRNA splicing factor 
SC35 in nuclear speckles (Figure 2C and Supplementary 
information, Figure S3F). Interestingly, WTAP deple- 
tion decreased the accumulation of both METTL3 and 
METTL14 in nuclear speckles (Figure 2B, 2C and Sup- 
plementary information, Figure S3G), whereas knock- 
down of either METTL3 or METTL14 had no effect on 
the localization of WTAP (Figure 2D and Supplementary 
information, Figure S3G). Furthermore, METTL3 deple- 
tion affected the nuclear speckle localization of MET- 
TL14 and vice versa (Figure 2B and 2C). Consistent 
with our observations on m6A methyltransferase activity, 
WT1 depletion did not affect METTL3 or WTAP local- 
ization (Figure 2B and 2D). Based on these observations, 
we conclude that WTAP is required for accumulation of 
METTL3 and METTL14 in nuclear speckles. 

Nuclear speckle-associating factors can be classified 
as RNase-sensitive or -insensitive factors based on the 
nature of the association with the nuclear speckles. To 
characterize the nature of the association of the WMM 
complex with nuclear speckles, we treated cells with 
RNase A. Similar to the mRNA processing factor SM 
(smith antigen), but unlike the splicing factor SC35, 
RNase A treatment resulted in a dramatic decrease of 
WMM complex accumulation in nuclear speckles, dem- 
onstrating that the WMM complex associates with RNA, 
which is essential for the recruitment and/or retention of 
the WMM complex in nuclear speckles (Supplementary 
information, Figure S3D-S3F). It is worth noting that 
METTL3 and METTL14 staining almost disappeared 
upon WTAP knockdown, but the overall protein levels 
of METTL3 and METTL14 did not change following 
WTAP depletion (Supplementary information, Figure 
S3B and S3C). It is likely that when WTAP is depleted, 
METTL3 and METTL14 no longer bind to RNA in the 
nuclear speckles, and become diffused within the entire 
cell, which reduces the staining density to levels below 
detection. In addition, unbound METTL3 and METTL14 
might easily get washed out during the permabilization 



and fixation procedures. A similar staining pattern in 
nuclear speckles has also previously been reported [15]. 

mRNA is the major RNA species bound by WTAP and 
METTL3 

We next performed Photoactivatable-Ribonucleoside- 
Enhanced Crosslinking and Immunoprecipitation (PAR- 
CLIP) [51] to identify the targets of WTAP and METTL3 
proteins. In brief, 4-thiouridine (4-SU), a photoreactive 
ribonucleoside analog, was incorportated into RNA tran- 
scripts, which leads to the thymidine-to-cytidine muta- 
tion in crosslinked sequences. Therefore, one unique 
feature of cDNA libraries prepared by PAR-CLIP is that 
the precise position of crosslinking can be identified by 
mutations residing in the sequenced cDNA. PARalyzer 
[36, 37] software was used to identify the preferred re- 
gions bound by WTAP and METTL3. 4 986 METTL3- 
and 922 WTAP -binding clusters were enriched and could 
be assigned into four RNA categories: mRNA, miRNA, 
large intergenic non-coding RNA (lincRNA) and pseu- 
dogene (Figure 3A). About 95% and 92% of the RNA 
bound by METTL3 and WTAP, respectively, represented 
mRNA, demonstrating that mRNA is the major substrate 
of the WMM complex. We further mapped the mRNA 
clusters into four non-overlapping transcript regions, 5'- 
UTR, CDS, 3'-UTR and intron, according to the human 
reference transcript annotation. After normalization of 
the number of clusters to the overall length of the dif- 
ferent transcript regions, the majority of binding sites of 
WTAP and METTL3 were identified in CDS and UTR 
regions (Figure 3B). 

To define the in vivo RNA recognition elements for 
WTAP and METTL3, the binding cluster data were ana- 
lyzed by HOMER, a suite of tools for motif discovery 
and next-generation sequencing analysis [38]. In this 
analysis procedure, the WTAP and METTL3 clusters 
were set as the target sequences, and a set of background 
clusters was generated with the BEDTools' shuffleBed 
program [39] to randomly shuffle regions of the same 
size as the clusters throughout the gene regions, with 
the parameter for HOMER motif length from 5 to 8. 
The motifs AGGACU (P = le-14) and UGUGGACU 
(P = le-13) were enriched in WTAP- and METTL3- 
binding clusters, respectively (Figure 3C). When we only 
included genes found in both WTAP- and METTL3- 
binding clusters, the highest scoring motif was UUAG- 
GACU (P = le-19) (Figure 3C). In addition, AGGAC 
(P = le-12), UGGAC (P = le-12), and AGACUAA (P 
= le-10) were also highly enriched in WTAP, METTL3 
and WTAP/METTL3 overlay clusters, respectively (Fig- 
ure 3D and Supplementary information, Figure S4B). 
This is in accordance with the reported consensus m6A 



www.cell-research.com | Cell Research 



<@ WTAP as m6A methyltransferase regulatory subunit 
182 



WTAP 



METTL3 



92.02 




0.90% 94.82%, 
5.06% 
2.02% 




2.17% 
2.57% 
0.44% 



WTAP 

UQUQGACU 

mam 



METTL3 



METTL3/WTAP 
OVERLAY 



I mRNA ■ pseudogene lincRNA BmiRNA 



B 



43.16% 



WTAP 

3.93% 



METTL3 




2.29% 



49.29% 



32.94% 



74% 



'35.17% 

H5.4I 

3'UTR ■ 5'UTR ■ CDS ■ intron 




METTL3 



WTAP 




METTL3 




WTAP 



METTL3 

Ref-gene- 
isoform 



Overlay with m6A 

5.45% 94.55% 
5.74% 94.26% 

61 .56% 38.44% 

i Single-isoform gene 
Multi-isoform gene 



siCTRL 
siWTAP 



IP:METTL3 




RNA:METTL3 




69.23% 



30.77% 



WTAP METTL3 overlay 
■ conservative ■ skipped 



Cell Research | Vol 24 No 2 | February 2014 



Xiao-Li Pinged a/. @ 



motif RRACH (R = G or A; H = A, C or U) [9, 1 0]. The 
high degree of similarity of the mRNA binding motifs 
between WTAP and METTL3 are consistent with the re- 
sults that WTAP and METTL3 form a complex. 

We next calculated the distance between WTAP/ 
METTL3 clusters and the m6A sites. To obtain the best 
coverage of m6A sites, we downloaded human m6A data 
released by Dominissini et al. [9] and Meyer et al. [10] 
and aligned the data to the hgl9 genome by Tophat2 [40]. 
We then used MACS [41] and Fisher test [9, 10] meth- 
ods to determine m6A peaks. The m6A peaks identified 
in this analysis were combined with the previously pub- 
lished m6A sites, and the duplicates were removed. Then 
we used the BEDTools' closestBed [39] to calculate the 
distance between WTAP/METTL3 clusters and m6A 
sites. We generated control clusters with BEDTools' 
shuffleBed [39] to randomly shuffle regions of the same 
size as the clusters. The result demonstrated a significant 
spatial correlation (P = 2.2e-16) with 74% WTAP- and 
69% METTL3 -binding clusters in CDS and UTR regions 
overlapping with m6A sites (Figure 3E and Supplemen- 
tary information, Figure S4C and S4D). These results 
strongly suggest a significant role of both WTAP and 
METTL3 in m6A modification. 

The results that WTAP is required for nuclear speckle 
localization of METTL3 and METTL14 and for m6A 
modification prompted us to test whether the absence of 
WTAP might influence the binding of METTL3 to RNA. 
Interestingly, we found a significantly reduced amount 
of RNA associated with METTL3 upon WTAP depletion 
(Figure 3F and Supplementary information, Figure S5B), 
suggesting that WTAP is likely responsible for recruit- 
ing the m6A methyltransferase complex to target RNAs. 



Taken together, these findings support the existence of 
a WTAP -METTL3 -containing complex that methylates 
RNA at specific sites in mammals. 

The WMM complex regulates transcription and 
alternative splicing 

Gene Ontology (GO) analyses of mRNA captured by 
METTL3 and WTAP using PAR-CLIP revealed that a 
large proportion of mRNA species associated with both 
proteins were derived from genes involved in transcrip- 
tion and RNA metabolism/splicing (Supplementary 
information, Figure S5A). To further investigate the po- 
tential roles of the WMM complex in RNA metabolism, 
transcriptome analyses were performed to compare the 
genome-wide gene expression and alternative splicing 
events between wild-type cells and WTAP- or METTL3- 
deficient cells. The transcriptome analyses revealed a sig- 
nificant overlap in differentially expressed genes (DEGs) 
between WTAP- and METTL3 -deficient cells (Supple- 
mentary information, Figure S5C). Consistent with the 
GO analysis of the PAR-CLIP data, GO analysis of the 
overlapping DEGs from WTAP- and METTL3 -deficient 
cells demonstrated that most of these genes were related 
to transcription and RNA processing (Supplementary 
information, Figure S5D), further emphasizing the poten- 
tial significance of WTAP and METTL3 in RNA metabo- 
lism. 

Interestingly, both METTL3 and WTAP knockdowns 
resulted in a pronounced variation in isoform numbers 
with about half of the variations (2 296) present in both 
WTAP- and METTL3 -deficient cells (Supplementary 
information, Figure S5E). When comparing the 2 296 
genes with changes in isoform numbers with the 410 



Figure 3 METTL3 and WTAP bind m6A consensus sequence in mRNA and affect gene expression and alternative splicing. (A) 
Percentage of various RNA species bound by WTAP and METTL3 based on PAR-CLIP analyses. The WTAP- and METTL3- 
binding clusters were identified by PARalyzer algorithm. Annotation of clusters was based on the human Ensembl gtf file (ver- 
sion 72, hg19). The majority of binding sites of WTAP and METTL3 were located in mRNA. (B) Pie chart depicting the distri- 
bution of binding clusters in mRNA based on PAR-CLIP sequence clusters for WTAP and METTL3 after normalization of the 
overall length of the different transcript regions. Length of these different transcript regions was extracted from Ensembl an- 
notations and the distribution percentage of clusters in these regions were normalized by their length. (C) Enriched sequence 
motif analysis of binding clusters indentified by PAR-CLIP. Upper panel, WTAP-binding motif ( P = 1e-14); middle panel, 
METTL3-binding motif (P = 1e-13); lower panel, binding motif obtained when only genes found in both WTAP- and METTL3- 
binding clusters were included (P = 1e-19). Binding motifs were computed by the HOMER program. (D) Venn diagram of the 
overlapping genes with binding clusters of WTAP and METTL3 in the PAR-CLIP samples. (E) Percentage of WTAP/METTL3 
clusters in CDS and UTR regions overlapped with m6A sites. (F) HeLa cells were transfected with siCTRL or siWTAP and 
Myc-METTL3 for 48 h as indicated. The cell lysates were then subjected to PAR-CLIP using anti-Myc. The pulled down RNA 
products in the RNA-METTL3 complex were labeled by Biotin and detected by Biotin chemiluminescent nucleic acid kit. (G) 
Percentage of WTAP- (711 multi-isoform and 41 single-isoform) and METTLE3- (3 155 multi-isoform and 192 single-isoform) 
binding mRNAs derived from single-isoform or multi-isoform genes and the reference Ensembl genes of human (P = 2.2e- 
16, Fisher test). (H) Percentage of constitutively or alternatively spliced exons adjacent to intronic binding clusters of WTAP 
(left), METTL3 (middle), overlap of WTAP and METTL3 (right). Supportive data were included in Supplementary information, 
Figures S4 and S5. 
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genes bound by the two proteins (Figure 3D), 99 genes 
were overlapped (Supplementary information, Figure 
S5E). Furthermore, assigning genes with METTL3- or 
WTAP-binding clusters to either single- or multi-isoform 
genes revealed that the majority of mRNA species bound 
by WTAP or METT3 were derived from multi-isoform 
genes (Figure 3G; P = 2.2e-16, Fisher test). Refining this 
observation further, assignment of all intronic clusters 
to either constitutively or alternatively spliced exons re- 
vealed that these clusters are significantly closer to alter- 
natively spliced exons compared to constitutively spliced 
exons (Figure 3H, P = 2.2e-16, Fisher test), further in- 
dicating a role of the WMM complex and mRNA m6A 
modification in RNA splicing as proposed previously [9]. 

The WMM complex is required for tissue differentiation 
in zebrafish 

The important role of the WMM complex in the regu- 
lation of m6A modification, gene expression and mRNA 
processing in cultured cells prompted us to investigate 
its functions at the organismal level using zebrafish em- 
bryos as the model system. Human and zebrafish WTAP 
proteins share 8 1 % identity with a high degree of conser- 
vation in N-terminal region (Supplementary information, 
Figure S5F). Whole-mount in situ hybridization (WISH) 
showed that WTAP and METTLS were maternally ex- 
pressed from the 4-cell stage and ubiquitously expressed 
through early embryogenesis, with enriched expression 
in the brain region at 36 hpf (hours post fertilization) 
(Supplementary information, Figure S6A). To determine 
whether WTAP and METTL3 are functionally involved 
in early embryonic development, WTAP and METTL3 
were knocked down in zebrafish embryos using antisense 
morpho linos (MOs). Embryos injected with WTAP MO 
displayed multiple developmental defects at 24 hpf in- 
cluding smaller head and eyes, smaller brain ventricle, 
and curved notochord, while embryos injected with 
METTLS MO were only modestly affected (Figure 4A). 
However, simultaneous knockdown of these two genes 
caused much more severe defects in embryos compared 
to parallel control and embryos injected with individual 
MOs (Figure 4A). GFP expression of the constructs of 
GFP-WTAP and GFP-METTL3 was markedly decreased 
by co-injection of corresponding MO, verifying efficient 
knockdown (Supplementary information, Figure S6B). 
WISH showed that the expression of neuro-ectoderm 
marker goosecoid, and mesoderm marker no tail at the 
50% epiboly stage was not changed, suggesting normal 
germ layer formation in the morphants (Supplementary 
information, Figure S6C). However, the expression of so- 
mite marker myod was increased in WTAP-, METTL3-, 
and double knockdown embryos (Figure 4B), indicating 



that the WMM complex might play an inhibitory role in 
muscle development. Although the proportion of apop- 
totic cells detected by the TUNEL assay was slightly 
increased in individual morphants, double knockdown of 
WTAP and METTL3 led to a striking increase in apop- 
tosis (Figure 4C), consistent with observations in human 
cells (Supplementary information, Figure S6D-S6F). 
The apoptotic cell death observed in WTAP MO-treated 
zebrafish embryos was rescued by either full-length or an 
N-terminal region of zebrafish WTAP (Figure 4D). Con- 
sistent with their role in mammalian cells, both METTL3 
and WTAP are required for the m6A methyltransferase 
activity in zebrafish embryos (Supplementary informa- 
tion, Figure S6G). These data suggests that the m6A 
methyltransferase is required for normal tissue differen- 
tiation in zebrafish. 

Discussion 

METTL3 homologs have been found in various spe- 
cies including S. cerevisiae (IME4), Drosophila (Ime4) 
and Arabidopsis (MTA) [23-26, 42]. The two consensus 
methyltransferase motifs CM I and CM II are present in 
all homologs. Single amino acid substitutions in the puta- 
tive catalytic domain of Ime4 lead to loss of m6A modifi- 
cation in mRNA and severe sporulation defects [42, 43]. 
In Arabidopsis, MTA is mainly expressed in dividing 
tissues, particularly in reproductive organs, shoot meri- 
stems and emerging lateral roots. Inactivation of MTA 
results in m6A modification defects and developmental 
failure beyond the globular stage [26]. 

Despite the identification of a multi-subunit complex 
possessing m6A methyltransferase activity more than 
15 years ago, METTL3 has remained the only known 
component of the mammalian complex [22]. There- 
fore, unveiling proteins associated with METTL3 in 
the methyltransferase complex is of key importance for 
understanding how m6A modifications in RNA might 
impact and regulate cellular functions. Previously, At- 
FIP37 and Mum2 were found to associate with plant and 
yeast METTL3 homologs, respectively [26, 28]. Like 
MTA, disruption of AtFIP37 also results in developmen- 
tal arrest at the globular stage [26, 44]. In the present 
study, we have identified WTAP, the human homolog of 
AtFIP37/Mum2, and a novel protein METTL14 as two 
new components of the human m6A methyltransferase 
complex. Like METTL3, METTL14 contains a catalytic 
domain and harbors catalytic activity [30]; while our data 
suggest that WTAP likely have very important regulatory 
functions, including tethering the m6A methyltransfer- 
ase complex to RNA and facilitating its accumulation 
in nuclear speckles (Figures 2B, 2C and 3F). The fact 
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Figure 4 The WMM complex plays essential roles during zebrafish embryogenesis. (A) Embryos injected with individual 
(WTAP-, METTL3-) or combined (METTL3+WTAP) MOs. The double knockdown showed more severe morphological defects, 
compared to other groups, at 24 hpf. Red arrows mark head and eyes, while blue arrows mark brain ventricle. The curve of 
the notochord is labeled by the double dashed lines. (B) Expression of somite marker myod in the morphants at 24 hpf. myod 
expression was increased in WTAP-, METTL3-, or double morphants. (C) Increased apoptosis (TUNEL assay) was observed 
in embryos injected with individual or combined MOs. (D) Overexpression of full-length or N-terminal but not C-terminal ze- 
brafish WTAP mRNA prevented apoptosis in zebrafish WTAP-UO embryos. Note that there is auto-fluorescence on the yolk. 
Anterior to the left and dorsal is up. Scale bar, 250 |im. Supportive data were included in Supplementary information, Figures 
S5F and S6. 



that m6A demethylases also localize in nuclear speck- 
les strongly suggests that a dynamic regulation of m6A 
levels in mRNA may occur in nuclear speckles [3, 14]. 
Furthermore, our studies in zebrafish revealed that MET- 
TL3 mRNA expression level changed during embryonic 
development (also observed in vitro in bovine embryos 
[45]), and that depletion of WTAP and METTL3 com- 
promised tissue differentiation (Figure 4 and Supplemen- 
tary information, Figure S6), strongly suggesting that 
m6A may play a key role in regulating organismal devel- 
opment. 

On average, there is 1 m6A modification per 2 000 



ribonucleotides [2]. Motif analyses of m6A peaks en- 
riched in both human and mouse cells have revealed that 
m6A modification occurs in highly conserved regions 
with a consensus sequence RRACH (R = G or A; H= A, 
C or U) [9-13]. The m6A modification occurs in intra- 
genic regions particularly in the 3'-end of CDS and the 
first quarter of the 3'-UTR [9, 10]. Our PAR-CLIP data 
demonstrate that the conserved WTAP- and METTL3- 
binding motifs both contain the consensus m6A sequence 
(Figure 3C), and their binding clusters in CDS and UTR 
regions significantly overlap with m6A sites (Figure 
3E), strongly suggesting a direct connection between the 
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WMM complex and m6A modification. 

Interestingly, without normalization of the overall 
length of the various transcript regions, the majority of 
WTAP- or METTL3 -binding clusters were located in 
introns (Supplementary information, Figure S4A), which 
suggests that WMM-dependent m6A modification may 
be coupled to gene transcription as previously proposed 
[1-3]. For the relationship of m6A and mRNA splicing, 
Dominissini and colleagues calculated the average num- 
ber of m6A peaks in multi-isoform and single-isoform 
genes, and assigned all m6A peaks to either constitu- 
tively or alternatively spliced exons based on the list of 
splicing events of all coding Ensembl genes. Their study 
revealed that m6A modification is significantly over-rep- 
resented in multi-isoform genes and alternatively spliced 
exons, proposing that m6A may affect RNA splicing [9]. 
Consistent with this observation, we found an enrichment 
of WTAP/METTL3 -binding clusters in genes with multi- 
isoforms and alternatively spliced exon junctions based 
on our PAR-CLIP+RNA-seq dataset (Figure 3G and 3H). 
Taken together, these results suggest that WTAP/MET- 
TL3 and thereby the m6A methyltransferase may influ- 
ence RNA splicing. However, a thorough understanding 
of the relationship between m6A and transcription and/or 
splicing awaits future investigations. 

Drosophilia female-lethal (2)D (FL (2)D, WTAP 
homolog) also interacts with the female-specific RNA- 
binding protein sex-lethal (SXL) [46, 47], VIR [46], Snf, 
U2AF, U2A50 and U2AF38 [47]. This together with the 
results showing that there is a difference in the numbers 
of genes containing METTL3- or WTAP-binding clus- 
ters might suggest that METTL3 is the core component 
of the m6A methyltransferase, whereas WTAP serves as 
its regulatory subunit. Based on the findings that WTAP 
interacts directly with METTL3 (Figure 1) and medi- 
ates nuclear speckle localization of both METTL3 and 
METTL14, we propose a model, in which METTL3 and 
METTL14 are anchored to mRNA by WTAP, allow- 
ing subsequent methylation of target adenosine residues 
(Figure 5). 

In summary, our findings identify WTAP as a regula- 
tory subunit required for formation of a functional m6A 
methyltransferase complex including METTL3 and 
METTL14, which play an important role in the regula- 
tion of gene expression and alternative splicing. Thus, 
the discovery and characterization of the RNA m6A 
methyltransferase complex, together with the recent ad- 
vances in the understanding of the demethylases FTO 
and ALKBH5, is an important step towards a more thor- 
ough understanding of the biological significance of the 
epitranscriptomic marker m6A. 




Figure 5 A schematic model of the WMM complex function. 
The WMM complex mediates methylation of internal adenosine 
residues in eukaryotic mRNA, forming N6-methyladenosine. 
WTAP binds to the m6A consensus RRACH motif of mRNA and 
recruits catalytic subunits METTL3 and METTL14. Then the 
METTL3-METTL14 complex carries out m6A methyltransferase 
activity in the m6A motif. 



Materials and Methods 

Cell culture, plasmids and antibodies 

Human cells were cultured in standard cell culture Dulbecco's 
modified Eagle's medium at 37 °C in a humidified incubator with 
5% C0 2 (v/v).The human METTL3 gene was cloned into pCMV- 
Myc (Invitrogen), S-protein/FLAG/SBP (streptavidin-binding 
protein) triple-tagged destination vector [48] or pGEX-5x-2 (GE 
healthcare). The human WTAP gene was cloned into p3XFLAG- 
CMV-14 (Sigma) or pProEX-HTb (Invitrogen). The human MET- 
TL14 gene was cloned into pcDNA3-HA (Invitrogen). The follow- 
ing antibodies were used: rabbit-anti-METTL3 (Abeam), mouse- 
anti-METTL3 (Abnova), mouse-anti-WTAP (Santa Cruz), rabbit- 
anti-METTL14 (Atlas), rabbit-anti-m6A (Synaptic Systems) and 
other antibodies included in Supplementary information, Data SI. 

Tandem affinity purification and mass spectrometry 

293T cells transiently transfected with triply-tagged METTTL3 
plasmid (SFB-METTL3) were used for tandem affinity purifica- 
tion. Cell lysates were incubated with streptavidin-Sepharose 
beads (Amersham Biosciences), and eluted with 1 mg/ml bio- 
tin (Sigma). The eluates were further incubated with S-protein- 
agarose (Novagen). The proteins bound to S-protein-agarose beads 
were subjected to SDS-PAGE and visualized by coomassie blue 
staining. The identities of eluted proteins were revealed by mass 
spectrometry analysis. The Q Exactive mass spectrometry data 
(Thermo Fisher Scientific) were searched against SwissProt hu- 
man database using 15 ppm peptide mass tolerance and 20mmu 
fragment mass tolerance. 

Immunoprecipitation 

293T cells transfected with the indicated siRNAs and/or DNA 
constructs were lysed in buffer (100 mM NaCl, 20 mM Tris- 
HC1 (pH 7.4), 0.5% NP-40, 1 mM PMSF, 1 mM Na 3 V0 4 , 1 mM 
P-glycerophosphate, 1 mM NaF and l x Cocktail), and subjected 
to immunoprecipitation (IP) followed by immunoblotting with the 
indicated antibodies. 
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The following antibodies were used in the study: mouse-anti- 
Flag (Sigma), rabbit-anti-Myc (Abeam), rabbit-anti-HA(Clontech), 
rabbit-anti-WTAP (Atlas), Anti-Rabbit IgG-HRP (Dakocytoma- 
tion), Anti-Mouse IgG-HRP (Dakocytomation). 

GST pull-down assay 

The human METTL3 gene was subcloned into pGEX-5X-2 
expression plasmid with GST-tag and the human WTAP gene 
was subcloned into pProEX-HTb expression plasmid with His 
tag. Recombinant GST-METTL3 and His-WTAP proteins were 
induced to express in E. coli strain BL21(DE3) and purified by 
FPLC using Bio-Scale Mini Profinity GST cartridge (Biorad) and 
Bio-Scale Mini Profinity IMAC cartridge (Biorad) as described 
in the commercial instruction manual. The quality of proteins was 
tested by western blot and coomassie brilliant blue staining. His- 
WTAP proteins were incubated with glutathione sepharose (GE 
Healthcare) to be pre-cleared in NETN buffer (100 mM NaCl, 1 
mM EDTA, 20 mM Tris-HCl pH 7.4, 0.5% NP-40). Then pre- 
cleared His-WTAP was mixed with GST or GST-METTL3 protein 
and incubated overnight at 4 °C with equal amount of glutathione 
sepharose beads. After washing the beads five times with NETN 
buffer, proteins bound to the beads were analyzed by 8% SDS- 
PAGE followed by immunoblotting with HRP-conjugated anti-His 
and anti-GST antibodies (Abeam). 

Immunofluorescence staining 

HeLa cells grown on coverslips were transfected with siRNA 
(50 nM) or DNA constructs as indicated for about 4 h, and 24 h 
later fixed with 4% paraformaldehyde on ice followed by permea- 
bilization with 0.1% Triton X-100 and 0.05% NP-40 on ice. After 
pre-blocking with 5% non-fat dried milk in TBST, coverslips were 
first incubated with primary antibody for 1 h and then with fluo- 
rescent dye-conjugated secondary antibody for 0.5 h and mounted 
with DAPI-containing mounting medium (Vector Laboratories, 
Burlingame, CA, USA) and visualized under immunofluorescence 
confocal microscope Leica TCS SP5 (Leica). The method used 
for immunofluorescence staining following RNase A treatment 
was adapted from Mayer et al. [49] and Sytnikova et al. [50]. In 
brief, cells grown on coverslips were permeabilized with 0.05% 
Triton X-100 in 20 mM Tris-HCl (pH 7.4), 5 mM MgCl 2 , 0.5 mM 
EDTA, and 25% glycerol; washed twice with PBS; and treated 
with RNase A (1 mg/ml) (diluted in PBS) in 37 °C for 7 min. After 
washing, cells were fixed with cold methanol (15 min) and stained 
with the indicated antibodies following the blocking procedure. 
The immunofluorescence images were digitally recorded by using 
confocal microscope. The following antibodies were used in im- 
munofluorescence staining: mouse-anti-Flag (Sigma), rabbit-anti- 
Myc (Abeam), rabbit-anti-HA (Clontech), mouse-anti-METTL3 
(Abnova), rabbit-anti-METTL14 (Atlas), mouse-anti-WTAP (Santa 
Cruz), mouse-anti-SC35 (Sigma), anti-mouse IgG-Cy3 (Sigma), 
anti-rabbit IgG-Cy3 (Sigma), anti-mouse IgG-FITC (Sigma), anti- 
rabbit IgG-FITC (Sigma). 

PAR-CLIP and RNA-Seq 

293T cells expressing Myc-METTL3 or Flag- WTAP were cul- 
tured in medium supplemented with 4-SU (Sigma), irradiated with 
365 nm UV light to induce crosslinking as described previously 
[51]. Immunoprecipitated protein-RNA complexes were separated 
by SDS-PAGE and then transferred to PVDF membrane; RNA- 
protein complexes were cut out from the membrane corresponding 



to the size of Myc-METTL3 or Flag- WTAP. Purified RNAs from 
RNA-Protein complex were subjected to the small RNA library 
construction and then sequenced with HiSeq 2000 system (Illumina 
Inc.). Bowtie 0.12.7 software [52] was applied to map sequencing 
reads against human (hgl9) genome with up to two mismatches 
allowed. To multi-mapped reads, mapped locations were only 
reported for those with the minimum number of observed mis- 
matches. PARalyzer software [36] was used to define METTL3- 
and WTAP-binding sites. Sequence motifs enriched in clusters 
were identified by HOMER, a suite of tools for motif discovery 
and next-generation sequencing analysis [38]. 

Total RNA was extracted from 293T cells transfected with 
siRNA using RNAzol (Molecular Research Center, Inc., USA). 
cDNA library was constructed using TruSeq RNA Sample Prep 
Kit and then sequenced with HiSeq 2000 system (Illumina Inc.). 
RNA-seq reads were aligned against the Ensembl genome (hgl9) 
using TopHat2 [40, 53]. The number of reads mapped to each of 
the Ensembl genes (version 72) was counted using the HTSeq 
python package [54]. DEGs were determined using the R-package 
DEGseq with the method MARS (MA-plot-based method with 
random sampling model), P value cutoff = 0.001 [55]. GO analy- 
ses of DEGs and PAR-CLIP data were performed using DAVID 
[56, 57]. Enrichment map of DEGs was constructed by Cytoscape 
2.8.3 software [57]. The FPKM (Fragments Per Kilobase of exon 
per Million mapped reads) value for each transcript was calculated 
by using Cufflinks version 2.0.0 and the differentially expressed 
transcripts were identified by Cuffdiff program in the Cufflinks 
software package [58]. Scripture, a method for transcriptome re- 
construction that relies solely on RNA-Seq reads and an assembled 
genome to build the transcriptome [59], was applied to calculate 
the isoform number for each gene. Variation in the isoform number 
between METTL3/WTAP-depleted and control samples indicated 
an alternative splicing event. 

Furthermore, a Bowtie library of non-redundant set of se- 
quences consisting of all possible junctions between the exons of 
each Refseq gene of human was generated and all the RNA-seq 
reads were aligned against this library with Bowtie 0.12.7 software 
[52], To the exon-exon junctions (EEJs) with read mapping, ac- 
cording to their continuous or discrete exon number, they can be 
divided into constitutively or alternatively spliced exons. All the 
WTAP/METTL3 clusters were assigned to these exons. 

MO-mediated gene knockdown in zebrafish embryos 

The antisense MO oligonucleotides targeting individual mem- 
bers of the WMM complex were synthesized (Gene Tools, LLC). 
The MOs (2.5-8 ng) were injected into embryos from one-cell 
to four-cell stage. Zebrafish WTAP and METTL3 were amplified 
from zebrafish cDNA library by PCR and sequentially subcloned 
into the pGEM-T vector and pEGFP-Nl (Clontech) to generate 
pWTAP-GFP or pMETTL3-GFP constructs for checking MO ef- 
ficiency 

Statistical analysis 

Two-tailed f-tests in grouped analyses of Prism5 software 
(GraphPad Software Inc., USA) were applied. P-values < 0.05 
were considered significant. 

Accession numbers 

Mass spectrometry data have been uploaded to Peptide Atlas 
under accession number PASS00296. PAR-CLIP and RNA-Seq 
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data have been uploaded to GEO database and can be accessed via 
accession number GSE50583. 
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